The Yangian symmetry and "{V8}" molecule 



Xu-Biao Peng*, Cheng-Ming Bai, Mo-Lin Ge 

Theoretical Physics Division, Chern Institute of Mathematics ,Nankai University, Tianjin 300071, 

P.R.China 



Abstract 

A new symmetry operator Q is introduced to re-describe the effective Hamilto- 
nian // 3 of the Heisenberg spin triangle in the {V6} molecule. The operator Q is a 
special form of Dzyaloshiky-Moriya (DM) interaction for the triangle with three 
spins itself. By extending Q to a system with four spins we present a new model 
comprised of four 1/2-spins and show that the ground state of the Hamiltonian // 4 
for the model may be with S = 1 by choosing some suitable exchange constants in 
H4. In analogy to the molecular {V6} where the interaction between two triangles 
in {V6} molecule is assumed also the DM type, we then give the theoretical mag- 
netization of the molecular model comprised of two four-spin systems. Supposed 
a single parallelogram model with 5 = 1 exists, the local spin moments are shown 
to be jqHb, JqI^b, Jo^b, jqVb respectively for the spin figuration. 
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1. Introduction 

The single molecular magnets(SMMs)have attracted much attention both for 
its scientific importance of studying fundamental issues and for its potential appli- 
cations. The magnetic molecules {V6} and {V15} provide us a good platform for 
exploring these issues, where the total spin is not large. There have been beautiful 



investigations on these molecules. iQlQ, SlH] 
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Figure 1: (color online)The structure model of {V6} molecule. The blue lines represents the two 
equal exchange interaction constants called J a , and the red line is the third exchange interaction 
constant J c . There is DM interaction H mter = A(Sa X SaO.v between the two triangles. 

For latter use, we will briefly introduce the structures of the molecular {V6} 
first. As was mentioned in the Ref.[Q], the molecular {V6} is the abbreviation of 
the magnetic molecule (CN 3 H 6 )4Na2[H4V 6 8 (P04)4{(OCH 2 )3CCH 2 OH} 2 ]- 14H 2 
Its structure is shown in FigQ] We see that the molecule consists of two pieces, 
each of which is an isosceles triangle. In each triangle, two of the 2-spin exchange 
constants are equal(shown in blue J a ~ 65K) and an order of the magnitude of the 
third(shown in red J c ~ IK). [3] And the experiment has shown that there is some 
special interaction called Dzyaloshiky-Moriya (DM) interaction[5] between the 
two pieces, whose Hamiltonian can be written as H inter = A(S A x S /i ') r [2] The 
operators S A and S A > are the total spin operators of two triangles, respectively, and 
the energy gap A is tiny. However, such an interaction can make contributions 
to the Landau-Zener-Stiickelberg (LZS) transition^ when the magnetic field is 
absent, which can be detected by measuring the magnetization of the molecular 
{V6}.i] 

In Ref.JlL 0, EL 0], 

the wave functions of molecular {V6} and the experimental 
measurements of the system have been shown clearly. However, from the theoret- 
ical point of view, each triangle of the molecular {V6} is formed by three spins, and 
the symmetry properties of the triangle desire to be investigated. In this article, 
we first focus on one piece of the molecular {V6}, and introduce a new symmetry 
operator Q to re-describe the triangular piece. We can see that the commutativity 
between such a new symmetry operator Q and the Hamiltonian H 3 will constrain 
the parameters in the Hamiltonian H 3 . The symmetry operator Q looks natural 
for the triangle model in the molecular {V6}, and further we shall extend such a 
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new symmetry to a four-spin system. Using the extended symmetry operator Q 
in the four-spin system, we establish its Hamiltonian H 4 whose model is usually 
a parallelogram. By analyzing the Hamiltonian H 4 , we find that the ground state 
can be with total spin S = 1 in some special cases, which has not been considered 
before. Based on this assumption, we have made a prediction on the magnetiza- 
tion of such a theoretical molecule comprised of two parallelograms with the DM 
interaction between them in analogy to the molecular {V6}.[|2|] To test the parallel- 
ogram model itself, we propose another molecular model and give the prediction 
of its local spin moments configuration. 

This article is organized as follows: In Secj2lwe will investigate the single 
Heisenberg spin triangle model in {V6} and introduce the new quantum number 
operator Q to describe it. We shall show that the operator Q represents a new 
symmetry in the three-spin system. In SecfJl we will extend the symmetry to a 
system comprised of four particles with each spin 1/2 and establish its Hamilto- 
nian H 4 . We demonstrate that the ground state of H 4 may be with total spin 5 = 1 
in some special cases. Besides, we make up a theoretical molecular model called 
{V8} in which there is the DM interaction between the two parallelograms and 
the prediction of its magnetization in experiment is made. In SecJH we propose 
another molecular model which contains only one parallelogram and give the pre- 
diction about its local spin moments configuration. In the Appendix, we shall 
show details about how the symmetry operator determines the Hamiltonian. 

2. the new symmetry in {V6} 

As was mentioned in the Sec[T] the molecular {V6} is comprised of two trian- 
gles with a tiny DM interaction between them. In this section, we will concentrate 
on only one triangular piece in the molecular {V6}. It has been testified that the 
triangle model is actually isosceles by the experiment. ^\ But in this section we 
would like to introduce a new symmetry operator Q to describe the triangle sys- 
tem. 

The Hamiltonian of the usual Heisenberg spin triangle can be written as: 

#0 = Jl2&l " $2 + ^23^2 - S3 + /13S1 • S3, (1) 

where Si, S 2 and S 3 are the spins of each particle, and the relationship among the 
exchange interaction constants J l2 , J 13 and 7 2 3 is unknown. If a magnetic field is 
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applied along z axis, then the term corresponding to the Zeeman Energy written 
as Hzeeman = ^B(S i z + S 2z + S 3z ) should be included in Eq.CQ.The Zeeman term 
can split the energy with different eigenvalues of S z , where S = Si + S2 + S3 is the 
total spin operator of the triangle. However, only the quantum numbers S 2 and S z 
are not adequate to describe a system with three spins. It is easy to see that there 
are two eigenstates corresponding to the quantum number S = l/2,S z = 1/2. 

The new symmetry operator Q we will introduce is written as Q = Y 2 , where 
the operator Y is a special form of the DM interaction in the Heisenberg spin 
triangle written as: 



It can be verified that such a new operator satisfies the commutation rules as 
[Q, S 2 ] = and [Q, S J = 0. So we will use the operator Q to represent a certain 
symmetry property of the three-spin system just like what we have done in the 
Hydrogen Atom. [7] (In Mathematics Physics, the operator Y is in fact a special 
form of the Yangian operator, see AppendixjAl) The operator Q can be viewed as 
a collective quantum number that describe the history besides S 2 and S z . If we 
take the set {S 2 ,S z , Q} to be the complete operator set of the system, the comrau- 
tativity [Q, H] = must be satisfied. From such an equation, we can easily get the 
relationship for the exchange interaction constants: J u = Jn- Fortunately, this 
relationship J\ 2 = J23, has been shown to exist in molecular {V6}[|81]. And with 
this symmetry, the Hamiltonian in Eq.CQ) can be simplified as: 



We emphasize that in the triangle model, the eigenvectors of the operator Q are 
nondegenerate, and the Yangian symmetry operator Q can uniquely determine 
the Hamiltonian of the system. So we can use the complete set {S 2 ,S Z ,Q} to 
determine the states in the triangular piece in the {V6} model described in SecfTJ 
By directly diagonalizing the matrix Q in the usual Lie Algebraic representation 
(see Appendix®, we get the two states with total spin S = 1/2, S z = -1/2 as: 



Y = /(Si x S 2 + S 2 x S 3 + Si x S 3 ). 



(2) 



Ht, - Jn{S\ • S2 + S2 • S3) + /13S1 • S3, 



\4> a ) = -^(ITU> + IUT>-2|4U», 
V6 

\<p p ) = -4=am>-uiT». 



(3) 



(4) 
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which are corresponding to the eigenvalues ~, -| of Q respectively. It can be 
verified that these states are the eigenstates of the Hamiltonian // 3 , and the corre- 
sponding levels are E\ _ 1 = - Jn and , = -f /13. 

2'~2 2' 2 

Hence, the special DM interaction operator Y plays an important role in the 
Heisenberg spin triangle in {V6}. Its square form Q is a new quantum number 
operator called Yangian symmetry operator. Using the symmetry operator Q, we 
can easily determine the form of the Hamiltonian, and directly obtain the proper 
states of the system. It is reasonable to extend such a Yangian symmetry to a four- 
spin system, and we will give a discussion about such a system in the next section. 



3. an extended {V6}-type molecular model 

3.1. a four-spin system determined by the Yangian symmetry 

For a system comprised of four particles with each spin 1/2, the most general 
Hamiltonian can be written as: 

4 

j>i=\ 

and the special DM interaction operator for the four-spin system is defined as: 

4 

Y = i^](S i xS J ). (6) 

j>i=l 

If there is the Yangian symmetry in the four-spin system, we still take the square of 
the DM interaction operator Q = Y 2 as a quantum number operator to characterize 
the system. We must emphasize that the DM interaction operator as Eq.© is 
identified with the special form of the Yangian symmetry operator in this four-spin 
system.(see AppendixjC]) Similar to the three-spin system in the Secf2l we shall 
get the constraint for the parameters in Eq.© from the commutativity [Q, H] = 0. 
To calculate the commutation relationship, we need get the eigenstates of Q first, 
and then let the Hamiltonian H share the same eigenstates with Q. However, it 
should be noted that in the four-spin system the eigenvectors of the operator Q are 
partially degenerate in the states with 5 = 1. We need introduce a more viable 
6 to form the degenerate states of the Hamiltonian H 4 as Eq.(fT9l) shown in the 
AppendixjC} Through the careful calculation shown in AppendixJO we obtain 
the Hamiltonian of the four-spin system as Eq.(l2"TT). In a special case(# = 0), we 
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(a) 



Figure 2: (a) (color online) The model of the four-spin system with the Yangian symmetry, the red 
line represent the exchange constants a\2, and the blue represent a\j,. (b) (color online) The model 
of the extended {V6}-type molecule, the interaction between the two parallelograms is H lnter = 
A(Sa x SA')y> which is totally similar to the {V6} molecule. 



can get the relationship between the exchange interaction constants in Eq.© with 
the Q- symmetry: 

a 12 = ^34, ao = a24, 

1 

«i4 = - (an + 2ai 3 ) , 

So the Hamiltonian can be written as a much simpler form 

Ha = «i2(Si • S 2 + S 3 • S 4 ) + ai 3 (Si • S 3 + S 2 • S 4 ) 

+ ^( a i2 + 2ai 3 )Si • S 4 + -(5ai2 - 2ai 3 )S 2 • S 3 . 

(8) 

The model of such a Hamiltonian is actually a parallelogram model[9] as shown in 
Figf2fa), where there are only two independent exchange constants called a\%, ai 3 . 

The eigenvalues and eigenstates of the Hamiltonian are as follows: 

#402,m = («12 + 2 fl 13)02,m, 
1 1 1 
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H^ 2 hm -- 


1 5 2 

= (gfll2 - ^«13)<Ai, m , 


H^l m -- 


4 5 , 

= (~3«12 + g«13)<Al, OT 


H^ox) = 


= (-2a u + ^fli3)^o.o» 


H 4 i// 0Q = 


3 i 



where m is the quantum number of S z and the eigenstates are shown as Eq.(fT8l) 
in the AppendixjCl If a magnetic field is applied, the states with different quan- 
tum numbers m will split because of the Zeeman Effect. When the field is weak 
enough, we emphasize that the state with total spin 5 = is not always the ground 
state. In fact, the ground state depends on the values of the exchange constants. 
By choosing some suitable constants, we can get the ground states with total spin 
5 = 1. For example, when the constants satisfy the relationship a\% < 0,au > 0, 
and a 12 > -2a 13, the ground state is 

l^_!> = -L (3| Till) + 1 ITU) - I IITI) - 3| HIT)) ; 

V20 

When an > 0, a l3 < 0, and an < -2a\ 2 , the ground state is 

1^ _i> = -7=0 Till) - 3| ITU) + 3| IITI) - I HIT)) ■ 
V20 

In both cases, the total spin of the ground state can be S = 1 . 
Thus, we have got the Hamiltonian of the four-spin system determined by the 
Yangian symmetry. In a special case it is a model of parallelogram. It should be 
emphasized that in such a system, the ground state may be no longer with total 
spin S = 0. If we take some suitable constants, the ground state can be with 5 = 1. 
This is an interesting case never considered before, and we will focus on such a 
case and make some experimental predictions. 



3.2. the magnetization of the extended {V6}-type molecule 

In the molecular {V6}, we have shown that there is a Yangian symmetry char- 
acterized by the operator Q in each triangle. Meanwhile it has also been proved by 
the experiment|0] that there is a DM interaction between the two triangular pieces 
in {V6} that may cause the LZS transition. [6] The effect of the LZS transition can 
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be reflected on the magnetization of the molecular {V6}. 



It is natural to extend the LZS transition for {V6} to a molecular system com- 
prised of two parallelograms. The molecular model is shown in FigJ2b), which 
can be called {V8}. We also assume that there is a similar DM interaction written 
as H inter = A(S A x S A ')v between the two parallelograms. The operators S A and 
Sa> are the total spin operators of the two parallelograms, respectively. Obviously, 
it is complicated to calculate such a system for any spin. To simplify the situa- 
tion, we only consider the ground state of each parallelogram as in the molecular 
{V6}.[3] If the total spin of ground state is zero, it will make no sense to study the 
magnetization. However, if the total spin is one, we can make a prediction about 
its magnetization. The Hamiltonian of such a system will be written as: 

M> = hyB(t)(S Az + S A , Z ) + A(S A x S A ,) y , (10) 

where the spins S A = 1 and S A > = 1, and y is the gyromagnetic ratio. The 
magnetic field B(t) is along the z axis and varies with the time. The first term 
in Eq. (fT0l) is about the Zeeman Energy, and the second one is the special DM 
interaction leading to the LZS transition and is assumed to be tiny. By exactly 
diagonalizing the Hamiltonian ^ ', we can get the energy levels shown in Figf3j 
The DM interaction plays an important role only near the "crossing point" of the 
energy levels when the magnetic field is absent. Far away from the "crossing 
point", the system behaves as two independent particles with each spin 5 = 1 
and we just need to consider one particle to investigate the magnetization of the 
system. In the experiment, the measurement of the magnetization is related to 
the relaxation time of the magnetization, and here we assume the relaxation time 
is long compared with the experimental time so that the hysteresis effect can be 
observed. The general method to investigate the magnetization has been fully 
discussed in the Ref.[l], and we shall directly make use of the conclusion in the 
reference. Using the standard formula 

~fT = X W N'N(t)PN>N> ~ X Wnn ' WPNN 
dt N' N' 

given in the Ref.[l], where the meaning of p and W can be seen clearly later, we 
can get the magnetization dynamic equation of the spin 5 = 1 written as: 



dt \ poo / \ C 3 C 4 \ poo ) \ F 



d( ZnZ \ ( C x C 2 \( &\JE. (U) 
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Figure 3: The scheme of the energy levels of the system comprised of two parallelograms with a 
DM interaction between them, where A is the energy gap of the Landau-Zener tunneling. 



where M is the magnetization of the system, and p is the density operator. The 

magnetization is written as:M = -hy{S z ) = -hy(p ++ - p ) and the coeffiecients 

in Eq.dTTI) are as follows: 



Ci = 


--(WLo - W +0 ) - W. + 


-w + ., 


c 2 = 


i(W +0 - W- ) + W 0+ - 


Wo- + W. 


c 3 = 


^(W +0 - W-o), 




c 4 = 


~(W +0 + W. + 2W 0+ 


+ 2Wb_), 


E = 


X -(W-v - W +Q ) + W. + - 


W + -, 


F = 


^(W +0 + W-o). 





The W NN > can be taken as the transition probability from the state \N) to \N'}, 
where N,N' is the magnetic quantum number characterized as +, 0, -. If we take 
the one phonon process approximation, the transition probability can be replaced 
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byrliflU 



W N n> = 



A8 3 



where 5 = E N -E N > = KyB(f)(N - N') is the energy difference between two differ- 
ent levels. If a pulsed magnetic field is applied, then a hysteresis loop(see Fig® 
can be obtained from the Eq. ([TT)) . where the magnetization has been normalized 
with M max = hy, and all the parameters in the Eq. (fTT|) are taken as the values of 
the {V6} molecule provided in the Ref.|2|]. We see that it is the usual hysteresis 
loop of a molecular magnet with the spin 5 = 1. 



However, in the {V8} model when the magnetic field varies to the vicinity of 
the "crossing point", the LZS transition must be considered. The Hamiltonian is 
written as Eq.dlOl). and the off-diagonal elements can not be ignored any more. By 
diagonalizing the Hamiltonian, the exact energy levels are: 



£2,0 
E2+2 



E\fl - £0,0 - 0; 

+ V# 2 (0 + A 2 ; 
1 



Eim = ±-= 



V2 
1 

V2 



5iu 2 B 2 (t) + 3A 2 + VV# 4 (0 + 30// 2 5 2 (0 + A 4 



5/u 2 B 2 (t) + 3A 2 - V 9n A B\t) + 30// 2 5 2 (0 + A 4 



where the subscript (/, m) is corresponding to the state \l, m) in the strong magnetic 
field in the -z direction. Here, I is the quantum number of total spin, and m is the 
magnetic quantum number. For simplicity, we just choose the nearest three levels 
E2 ,-1, E 2 ,o, £2,1 in a weak field to investigate the effect of the off-diagonal elements. 
These three levels can be viewed as the eigenvalues of the Hamiltonian with the 
matrix form 



( B(t) 

A 
V2 





A 
V2 



A 
V2 





A 
V2 

-B(t) 



(12) 



whose eigenstates are 

\E + ,t) 



1+cosyS sin/? l-cos/3 

2 '^r 2 
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Figure 4: (color online) The magnetization of the four-spin system in a low temperature. The blue 
solid line is the magnetization, and the red dashed line is the density matrix poo- The function of 
the magnetic field varies with the time as B{t) = 10 sin(f). 



\r a ,'- sin/3 sin/3 
\E ,t) = | — , cos/3, 



V2 



vr 



\E-,t) = 



l-cos/3 sin/3 l+cos/3 



2 V2 2 

where the angle /3 is determined by the formula cos /3 



B(f) 



(13) 

and the cor- 



responding eigenvalues are Eq = 0, E ± = + V^ 2 (t) + A 2 . Following the 
method shown in Ref.[l] and assuming the Landau-Zener tunneling adiabatic, we 
can get the magnetization for spin 5 = 1: 



M 



-hy{S z ) = -Hy(p ++ S + z + +pwS™ +p__5r) 



It is easy to calculate the matrix elements of S z ,and the Eq.(fT4l) can be simplified 
as M = hy cos fi-n. Here, n is defined as n = p__ -p++, and obeys the same Bloch 
equation as Eq. (fTT|) : 



d_l n{t) 
dt\ Poo 



Ci C 2 
C3 C4 



Poo / IF 



(14) 
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Figure 5: (color online) The magnetization of the system comprised of two four-spin system with 
DM interaction between them. The blue solid line is the magnetization, and the red dashed line is 
the density matrix poo- We can see the density matrix poo is not affected by the off-diagonal matrix 
of the Hamiltonian, and its magnetization has a similar behavior to the molecular IV61. lfTlll 

Combining the Eq.(fl4l) and (fl4l) . we can get the magnetization of the system as 
shown in FigJ51 where the applied magnetic field is taken as B(t) = 10 sin(f) vary- 
ing in a period t e [0, n]. Similar to the {V6} molecule, the Magnetization can be 
probed because of the LZS effect. 

In a summary, we have extended the properties of the molecular {V6} to a 
{V8}-like molecule comprised of two parallelograms. By investigating its magne- 
tization of the molecular model, we get the conclusion that the {V8}-like molecule 
has a similar behavior to the molecular {V6}. 



4. the local spin moment configuration of the four-spin system 

In this section, we will make up another molecule whose model can be treated 
as a parallelogram and give a prediction about its local spin moments configura- 
tion. 

To make up the theoretical molecular model, we need recall the model of 
{VI 5} which can be considered as an isosceles triangle in a low temperature first. 
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Figure 6: (a)(color online) The structure model of the molecular {V15}. The solid circles repre- 
sent the V 4+ = 1/2) ions. In a low temperature, only the triangle in the middle is considered. 
(b)(color online) The theoretical molecule whose model is a parallelogram satisfying the condition 
Eq.©. The upper and lower polygons are comprised of 2N particles so that their total spins in the 
ground state is zero in the low temperature respectively. The different colors represent the different 
coupling constants. 



The molecular {V15} is comprised of 15 V 4+ ions with each spin s = 1/2, and the 
ions are arranged in a quasispherical layered structure with a triangle sandwiched 
between two hexagons as is shown in the Figjfja). Each hexagon of {V15} consists 
of three pairs of strongly coupled spins with Ji ~ 800^.111211 Each spin of the V 4+ 
ions in the central triangle is coupled with the spins in both hexagons with J 2 = 
150K and 7 3 = 300^. 111211 resulting in a very weak exchange interaction between 



the spins within the central triangle with J = 2.44^.1130 When the temperature 
is low enough, the molecular {V15} has a much simpler approximation: the two 
hexagons can be omitted because each total spin is zero, and the only part we 
need to consider is just a simple Heisenberg spin triangle. [1^1 And it has been 
revealed that such a Heisenberg spin triangle model is an isosceles triangle with 
the relationship 7 13 > J n = Jz3 by measuring the local spin moment configuration 
in an NMR experiment. [4] So the Yangian symmetry also exists in the {V15} 
molecule in a low temperature, and it can be detected by measuring the local spin 
moments configuration. 

As the next step, we will extend such a molecular model to the four-spin 
system and give the prediction about the local spin moment configuration. The 
molecule we make up is shown in the Figjfjb), whose structure is totally similar to 
the {VI 5} molecule. The numbers of the ions in the upper and the lower polygons 
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are both even. And they are coupled to the special parallelogram in the middle, 
which results in some special exchange constants in the middle parallelogram sat- 
isfying the Eq.©. Similar to the {V15} molecule, such a molecule can be viewed 
as a single parallelogram described in Secf3l We can calculate the local spin mo- 
ments configuration of this model in the NMR experiment. The local spin moment 
is usually written as /i,- = -gn B {<f> g \S?\<f> g ), where the operator S ?(i = 1, 2, 3) is the 
z-component of the spin operator of the j-th particle, g = 2, and fi B is the Bohr 
magneton. The state \(f> g ) is the ground state of the system. Through direct calcula- 
tion, we can obtain the local moments of the four particles. It is highly nontrivial 
that the ground state of a four-spin system is with S = 1, so we mainly focus on 
such a case and predict its measurement in the experiment. When the exchange 
constants satisfy the relationship a 12 < 0, a i3 > 0, a 12 > - 2ai3, the ground state 
is _ v and the corresponding local spin moments are jq/^b, JqVb, Jq^b, JqMb, re- 
spectively. So the observations are jxfJ-B, j^/^b when the NMR experiment is per- 
formed. In the other case with an > 0, < 0, an < -2an, through similar 
calculation, we can find that the observations are also -^Hb, Jo^b- Hence, we can 
see that the measured local spin moments are always jq^b, JqHb when the ground 
state is with S = 1. So we can tell whether there is the Yangian symmetry by 
measuring its local spin moments. 

5. summary 

By investigating the experimental results on the molecular {V6}, we introduced 
a special form of the DM interaction operator to describe the symmetry of the 
molecule. Meanwhile we have extended the symmetry operator to the four-spin 
system with each spin 1/2. We find that the ground state may be no longer with 
total spin S = in a four-spin system determined by the Yangian symmetry. By 
choosing suitable interaction constants, the state with spin 5 = 1 can be the ground 
state. When the ground state is with spin S = 1, we extended the DM interaction 
in the {V6} molecule to our new model {V8}. By investigating its magnetization, 
we find that it has a similar behavior to the {V6} molecule. At last, we have 
proposed a theoretical molecular model to make a prediction about its local spin 
moments configuration which can be measured by the NMR experiment. 
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A. the mathematics of Yangian 



Yangian algebras were established by Drinfel'd.lla O, UM A Yangian is 



formed by a set {I, Y} obeying the commutation relations: 

[[Ya, [Yfi^v]] - Ua, [Y^, Y v ]] = h 2 a^ va py{I a ,Ii3,Iy}, 

U7a, Y,], \I<r, YrrJ] + [[^, Y x ], [h, Y,]] 

— h (tt AiivaftyC crxv & crxvafiyC AyLv)\Iai Ifli Yy\, (15) 

where the set of lx forms a simple Lie algebra characterized by c^ v and the re- 
peated indices mean summation. The definitions of a^ va js y and [xi, Xj, xu) were 
given in the Ref . [| 1 6T| . 

1_ 

Q-Ajivafly — "^CAao-CfifiTCvypCcrrpi 
{X\,X2,X^\ = XjXjXk. 

In the system comprised of n particles, if the operator I is taken as the total spin 
operator S = £" =1 S» the Yangian operator Y can be realized in terms of 

n n 

Y^^ + i^xSj), (16) 

where n is the number of the particles, and S, is the spin of the z-fh particle. It 
can be verified that the operators S and Y satisfy the Eq.(fT5l) and they form a 
Yangian. For a system with the Yangian symmetry, we introduce the symmetry 
operator Q = Y 2 to characterize it. So the commutation relation [Q, H] = must 
be satisfied, which means that both the operator Q and the Hamiltonian share the 
same eigenvectors. And the Hamiltonian of the n-spin system is usually written 
as: 

H = Y j a ij S r Sj, (17) 

«*</' 

where a, 7 is the exchange constant between the z-th particle and j-th particle. We 
note that the Hamiltonian is a symmetric matrix in the representation of Lie al- 
gebra. To satisfy the commutation relation, the symmetry operator Q is also de- 
manded to be symmetric in the representation. 
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B. The three-spin system with the Yangian symmetry 



When the system is a triangle system, the Yangian operator Y is taken as 
Eq.(fT6l). and the Hamiltonian is as Eq.(fT71). where n is taken as n = 3. We need 
to calculate the eigenvalues and eigenstates of Q first. The operator Q can be 
calculated as: 



Q = 2^ ul jS 2 j + 2(miM2Si • S2 + U2U3S2 • S3 + • S3) 



!=1 



+2i(ui - u 2 + u 3 )S 1 • (S 2 x S 3 ) - {S^Sj + SjSj + 
-(Si • S 2 + S 2 • S 3 + Si • S 3 ) + 2S 2 (S 2 • S 3 ) + ISjiSi ■ S 2 ) 
-25 ^ • S 3 ) - (Si • S 2 + S 2 • S 3 + Si • S3) 2 
+2[(S 1 -S 2 )(S 2 -S 3 ) + (S 2 -S 3 )(S 1 -S 2 )]} 

And the usual Lie algebraic bases are as follows: 



103 _4 > 



|03 _ 1 > 



h l) 

2'2 

h 3> 

2'2 



1 I 

2 • 2 



) = 



101,1) 



I0i -1) 



I III), 

-Ul Til) + 
V3 

-Ul TTI) + 
V3 

I TTT), 

-4=0 ITT) + 
V6 

4=d TID + 
V6 

-4=0 ITT) + 
V2 

4=(l Til) + 
V2 



iu)+iiim 

TIT) + 1 ITT)), 

TIT) -2| TTD), 
ITD-21 ||T», 
TIT)), 
ITD). 



It can be easily testified that 



3 1 

1 



2101 „,) = [7( M l+ M 2 + M 3) + XWlW2-"2"3-MlM3-7]l0'l m ) 
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V3 

2~("1 - "2 + 1)("3 + l)l^l, m >, 

Q\<f>i- >m ) = -^l-"2-l)("3-l)l0i^ + [^("l-"2) 2 + ^-^]]0^ m >. 

As mentioned in the AppendixjAl the Yangian operator must be symmetric, 
from which we can get the relationship w 2 = u\ + u?. If we take u\ = u 2 = 
ut, = 0, we can see that the Yangian operator is reduced to the special form of the 
DM interaction operator as Eq.©. In this special case, we can easily obtain the 
eigenvalues and the eigenvectors of Q as: 

Q\<f>l,J = -l|0 i>m >, 
where m is the magnetic quantum number, and the eigenvectors are 

= l<% m >, 

= -%^ + \\^y 

When m = -1/2 the states with total spin S = 1/2 are reduced to the Eq.© 
and Eq.©. We can see that the operator Q just mixes the states with the same 
numbers of S 2 ,S Z . To calculate the commutativity between the operator Q and 
the Hamiltonian, we demand the Hamiltonian share the same eigenstates as Q, 
i.e., 

H\ilA ) = E° |<A7 ), 



By direct calculation, we can get the relationship an = a 2 3 among the exchange 
constants in the Hamiltonian ([T71) . And the corresponding energy levels can be 
obtained as:£f , = \J n - J 12 , , = -|/ 13 , E h _s = J -f + J -f. 

2' 2 2' 2 2 2 
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C. The four-spin system with the Yangian symmetry 

For a system comprised of four spins, the similar process will be repeated. 
The Yangian operator Eq.(fT6l) and the Hamiltonian Eq.(fT71) are taken as those 
with n = 4. The operator Q in the four-spin system can be calculated as: 

4 4 

Q = J]u]Sj + 2 £ UiUjSrSi 

i=l i,j(i<j) 

+2i[(Ui - U2 + i<3)Si • (S2 X S3) + (U\ - U 2 + W/OSi • (S2 X S4) 

+(«1 - H3 + U 4 )Si ■ (S3 X S4) + (U 2 - «3 + M 4 )S2 ' (S3 X S4)] 
4 4 4 

ijXK/) ij(i<7) i,j(i<j) 

+2[S 2 l (S 2 • S 3 + S 2 • S 4 + S 3 • S 4 ) + 5^(S 3 • S 4 - Si • S 3 

—Si • S4) + S 3(81 • S2 — Si • S4 — S2 • S4) 

+5^(S 1 -S 2 + Si-S 3 + S 2 -S 3 )] 

+2[(S 2 • S 3 )(Si • S 2 ) + (Si • S 2 )(S 2 • S 3 )] 

+2[(S 2 • S 4 )(Si • S 2 ) + (Si • S 2 )(S 2 • S 4 )] 

+2[(S 3 • S 4 )(Si • S 3 ) + (Si • S3XS3 • S 4 )] 

+2[(S 3 • S4XS2 • S 3 ) + (S 2 • S 3 )(S 3 • S 4 )] 

+2[(Si • S 3 )(S 2 • S 4 ) - (Si • S 4 )(S 2 • S 3 ) + 3(Si • S 2 )(S 3 • S 4 )]} 

The usual Lie algebraic bases are as follows: 

102,-2) = HID, 

\<f> 2 ,-i) = j(imi> + iiru> + mu>+uim 

I02.o> = 4=(l TTU>+ tm>+ TUT>+ iUT>+ UTT>+ iTU». 
V6 

102.1) = ^(imt> + IUTT> + l4TTT>+Tm», 

102.2) = I TTTT), 
1 

2 



ill) = -HTm) + imT)-IUTT)-IITTT)), 



|0l 0) = -UlTTU>-IUTT», 
V2 



18 



I<_!> = ^(IUU) + UTU>-IUU)-UUT)-X 
I0?x> = 4=(IUH>-UTH», 

V2 

i0? >o ) = ^(i uu> + 1 rut)- 1 mr>-um», 
i0?_i> = -Uiuu>-um», 

i^!) = 4=(irm>-imr», 
V2 

I0?, o > = ^(iuu> + um>-imT>-iTum 

l^_i> = -UlUU>-UUT», 

l<o) = ^[2(1 TTU) + I UTT»- (I UU) 

+1 IUT> + ITUT> + Um»], 
I0o,o) = ^(1 UU) + UUT>- 1 TUT) -um». 

Then the action of Q turns out that 

3 2 1 2 5 

<2l02,m) = + "2 + "3 + "4) + ^("1 - U 2 ) - - 

1 2 1 2 

+ -(«! + i< 2 - "3 - "4) + ^("3 - "4) ]l02,m); 

1 9 1 

Gl01,m> = [^("l +W2 + M 3 +M 4)- j + 4^ Ml + "2 - "3 - "4) 2 ]|0^ m ) 

1 2 1 3 

=(m1 - M 2 + 1)(«3 + U4 + 2)|0! > + — (w 3 - w 4 + - "2 - 2)10! >; 

V2 V2 

m ) = f( m1 - "2 + 1)(«3 +U4 + 2)|0} m ) + [-(« 3 + i/ 4 ) 2 + t(«3 ~ U A f 

V2 '24 
3 1 

+ 4O1 - W2) 2 - l]l0 2 , m ) + ^("1 - "2 - 1)(«3 - "4 + l)l0l >m ); 

2101 m > = -^=(U 3 - U A + 1)(«1 - U 2 ~ 2)|0} m ) + i(«i - M 2 - 1)(«3 - "4 + 1)101 m ) 
V2 2 

113 

+[-Oi + i/ 2 ) 2 + ^Oi - w 2 ) 2 + ^(w 3 - w 4 ) 2 - i]l0i, m ); 
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2l0o,o> = {^("i + u 2 -u 3 -u 4 - 2)(«i + m 2 - "3 - u A + 2) + - u 2 ) 2 + (u 3 - u A f - 2]}|0 OO ) 

V3 2 

Z-( M 1 - "2 + 1)(«3 - "4 + l)l0 O ,o); 

V3 3 

QI^O.o) = 2~( M 1 - "2 + 1)(«3 - "4 + 1)I0O,()) + ^[("1 - M 2) 2 + ("3 ~ U A f - 2]|0^ o >; 

However, when the operator Q is required to be symmetric, we get the conclusion 
U\ = u 2 = w 3 = U4 = 0. So the first term of the Eq. (fT6l) must vanish, and the the 
Yangian operator is identified with special form of the DM interaction operator. 
In the representation of the Lie algebra, the eigenvalues and the eigenstates of the 
Q are: 



Q 



where the eigenstates are 



' 4>2,m ' 


f 


















v ^00 > 


V 



_5 

2 











_1 

2 













2 

























-1 





V 





-3 ) 



<p2,m 
Tl,m 

^00 
^00 J 



l02,-2> 
ltf.-l> 



l*Ll> 



^(-1 UU> + UTU> + IUU>-IU4T». 



1 



V20 
1 

V20 
1 



(31 UU> + I4TU>-IUU>-3| U4T», 



(I UU)-3| ITU) + 3| UTD-UUT)), 



:[(l TTU> + IUTT»-2(| UU>, 



l<Ao + ,o> 

l«Ao,o> = -l-(ITUi) + liiTT)) + (ITIiT) + UTTi»l- 



2V3 

+1 mr» + d tUT>+iim»L 
1 



(18) 



We can see that the eigenstates with S = lof<2 = -l/2 are partly degenerate. 
Here, we need another viable 8 to form the eigenstates of the Hamiltonian with 
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S = l,S z = -l,Q = -1/2 as: 



l^) =cos||a}_ 1 )-sinf|4_ l ) 

=sinf|or}_ 1 > + cosf|^_i> d9) 



With the commutativity between the operator Q and the Hamiltonian, we de- 
mand these six states |0 2 _ 2 >, l<A}'_i>> I^i,_i>. I<A?'-i>> I^o,q>> l<Ao,o> to be the eigen- 
states of the Hamiltonian, from which we can get the relationship among the in- 
teraction constants. The relationship is as follows: 

^24 = ^ ( ai2 + ^° 13 _ Q ^ ' 

a u = -(a n + 2a l3 ), 

«23 = ^ (2ai2 - 2a n + 3a 34 ) , (20) 

and the constants a l2 , 034 and a 13 satisfy the equation: 

2 (a n a 34 \ . 2 9 (5 5 \ 1.^/2 „ 8 \ n 
cos -y— 2~J + Sin 2\2 an ~ 2 a34 ! + 2 Sm \3 an ~ 3 au ! = 

Hence, the general four-spin Hamiltonian with the Yangian symmetry is written 
as: 

H = a 12 S 1 .S 2 + a3 4 S3-S 4 + a 13 S 1 .S3 + i(a 12+ 2a 1 3-a 34 )S 2 -S 4 



+\{a n + 2a 13 )Si • S 4 + ^(2a 12 - 2a 13 + 3a 34 )S 2 • S 3 , 



(21) 



where a n , an and a 3 4 satisfy the Eq.(l2TI). If we take the viable 6 to be zero, then 
the Eq.(l2"0"l)and(l2"TT) reduced to the Eq.© leading to a special parallelogram model, 
and the Hamiltonian Eq.(l2TI) reduces to Eq.©. 
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